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We map out the detuning-magnetization phase diagram for a "magnetized" (unequal number of 
atoms in two pairing hyperfine states) gas of fermionic atoms interacting via an s-wave Feshbach 
resonance (FR). The phase diagram is dominated by coexistence of a magnetized normal gas and 
a singlet paired superfluid with the latter exhibiting a BCS-Bose Einstein condensate crossover 
with reduced FR detuning. On the BCS side of strongly overlapping Cooper pairs, a sliver of finite- 
momentum paired Fulde-Ferrell-Larkin-Ovchinnikov magnetized phase intervenes between the phase 
separated and normal states. In contrast, for large negative detuning a uniform, polarized superfluid, 
that is a coherent mixture of singlet Bose-Einstein-condensed molecules and fully magnetized single- 
species Fermi-sea, is a stable ground state. 



Recent experimental realizations of paired superfluid- 
ity in trapped fermionic atoms interacting via a Fesh- 
bach resonance (FR) 0, 0] have opened a new chap- 
ter of many-body atomic physics. Almost exclusively, 
the focus has been on equal mixtures of two hyperfine 
states exhibiting pseudo-spin singlet superfluidity that 
can be tuned from the momentum-pairing BCS regime 
of strongly overlapping Cooper pairs (for large positive 
detuning) to the coordinate-space pairing Bose-Einstein 
condensate (BEC) regime of dilute molecules (for nega- 
tive detuning) 

In contrast, s-wave pairing for unequal numbers of 
atoms in the two pairing hyperfine states has received 
virtually no experimental attention and only some re- 
cent theoretical activity 0, 0, IE 0> 0, @ • Associating the 
two pairing hyperfine states with up (f) and down (J.) 
pseudo-spin a, the density difference Sn = n-\ — n.j is 
isomorphic to "magnetization" m = Sn and the corre- 
sponding chemical potential difference Sfi = — fi± to a 
purely Zeeman field h = Sfi/2. 

This subject dates back to the work of Fulde and Fer- 
rell (FF) 10] and Larkin and Ovchinnikov (LO) [H 
who proposed that, in the presence of a Zeeman field, 
an s-wave BCS superconductor is unstable to magne- 
tized pairing at a finite momentum Q ss fcp| — &fj. with 
kp a the Fermi wavevector of fermion a. This FFLO 
state, which remains elusive in condensed matter systems 
where it is obscured by orbital and disorder effects, spon- 
taneously breaks rotational and translational symmetry 
and emerges as a compromise between competing singlet 
pairing and Pauli paramagnetism. 

Thus atomic fermion gases (where the above delete- 
rious effects arc absent), tuned near an s-wave FR, are 
promising ideal systems for a realization of the FFLO 
and related finite-magnetization paired states, that can 
be studied throughout the full BCS-BEC crossover. 

In this Letter, we map out the detuning-magnetization 
phase diagram (FigQJ of such paired superfluids. We 
find that for positive detuning 8 and arbitrarily small to, 
the system phase-separates into a magnetized normal gas 
(N) and a singlet-paired BCS superfluid that exhibits a 
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FIG. 1: (Color Online) Detuning, 5 - population difference, 
m/n — (nj — ni)/{n^ + ny) phase diagram (for coupling 7 = 
0.1) in (a) displaying "normal" (N), magnetized superfluid 
(SFm), FFLO (thick red line) and SF-N coexistence states, 
(b) showing the FFLO wavevector Q(8) along the FFLO-N 
phase boundary, and (c) zoom-in on the FFLO state, stable 
only for 5 > <5* ~ 2.2eF. To the right of the dashed lines in 
(a) and (c), the SF-N coexistence undergoes a transition to 
SF-FFLO coexistence. 

BCS-BEC crossover with reduced 8. The FFLO state in- 
tervenes in a sliver on the boundary between this coexis- 
tence region and the N state. For large negative detuning 
a uniform magnetized superfluid (SFm), that is a coher- 
ent mixture of singlet Bose-condensed molecules and fully 
magnetized single-species Fermi-sea is a stable ground 
state. Our predictions of these states and transitions 
between them are testable via thermodynamics (qualita- 
tively modified by gapless atomic excitations inside the 
SFm and FFLO states), sound propagation (with zeroth 
sound velocity vanishing at the SFm-N transition), and 
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time of flight imaging (displaying density discontinuity 
and striking Bragg peaks associated with the finite mo- 
mentum pairing in the FFLO state). 

We now sketch the analysis that led to these results. 
A gas of fcrmionic atoms, a k a, resonantly interacting 
through a diatomic (closed-channel) molecule, b q , is de- 
scribed by a two-channel Hamiltonian [j| 



two of above pure states coexist as a mixture in fractions 
1 — x and x to be determined. 

The computation of the ground-state energy is simpli- 
fied by noting that Ec{h) — E G (0) — f£ m(h')dh', where 
Eq(0) is the well-studied fully-paired h = energy and 
m(h) = —dEa/dh is the atom species imbalance num- 
ber . We compute Eg by first neglecting the FFLO state 
(i.e., Q = 0), which, as we shall show, is only stable for 
a narrow window of parameters (see Fig.Q. Then, 
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where e k = k 2 /2m a , //f^ = /i±/i are the chemical poten- 
tials to impose atom number in hyperfine states |, I, or 
the total atom density n = + ri| + 2 ^2 q (b\bq) (imposed 

by fi, with n<j — J2 q airier ®ka)) an d density difference 
(magnetization) to — rij — (imposed by h). Here, So is 
the bare FR detuning, g is the FR coupling determining 
the resonance width and the system volume is unity. 

For a narrow FR (small <?), H can be accurately ana- 
lyzed by treating b q as a single-momentum [lfl Inl H2I 
c-number mode (bq) = bQS Qt Q with corrections small |l4| 
in powers of 7 = j(ep) = g 2 Af(ep) / 'cf, the ratio of the FR 
width to Fermi energy, with J\f(ep) = m 3 / 2 ^/ep j \Z2n 2 = 
Cy/eF the density of states at the Fermi energy ep = 
fc F /2m set by the total atom density n = gee 3 / 2 . To 
lowest order in 7, standard Bogoliubov analysis [15j gives 
the ground-state energy (h = 1): 
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E G = (H) = ( e -S- + 5 Q -2n)^-~ ]T> fe - e k ) 
+ ^^ kT e(-£ kT )+£ ki G(-^)], (2) 



where E^ a = Ek =F (h + k • Q/2to ) is the excitation 
spectrum for a hyperfine state <r, with "gap" Aq = gbQ 
and E k = (e^+A 2 ,) 1 / 2 , e k = fi+Q 2 /8m a , and Q(x) 

the Heaviside step function. The corresponding ground 
state is of the BCS form, but with pairing limited to 
momenta k satisfying E^ a > 0. 

The phase diagram is determined by minimizing Eg 
over Q and Aq at fixed average total density n, pop- 
ulation difference m, and physical detuning S — 60 — 
g 2 J2k^/^ e k (determined by the 2-body scattering am- 
plitude). The competing ground states are: (i) a normal 
Fermi gas (N) with Aq = 0, (ii) a non- "magnetic" fully 
paired BCS-BEC superffuid (SF) with A Q ^ 0, Q = 0, 
and to = 0, (iii) a magnetized partially paired superfluid 
(SFm) with Aq 7^ 0, Q = 0, and m^0, and (iv) a mag- 
netized, finite-momentum paired superfluid (FFLO) with 
Aq 7^ 0, Q 7^ 0, and m 7^ 0. Anticipating the existence 
of first-order transitions, across which m, n are discon- 
tinuous, in order to guarantee a solution everywhere it 
is essential to also include phase-separated states where 



m(h) = -c@(h - A) [{fi + y/h 2 - A 2 ) 3 / 2 



-(fi - y/h 2 - A 2 ) 3 / 2 6(/i - Vh 2 - A 2 )] . (3) 

For positive detuning S ^> 6F7 1 / 2 , appropriate in the BCS 
and throughout most of the crossover regimes, A -C /x 
and the density of states inside Eg(0) can be well ap- 
proximated by a constant jV(/i), giving: 
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m(h')dh'. 
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For small h -C /1 the species imbalance contribution to 
Eg is well approximated by J Q h m(h')dh' ~ N(n)Q(h — 
A)[hy/h 2 -A 2 - A 2 cosbT 1 (h/A)]. For < h < 
Abcs/2, Eq exhibits a single minimum at a standard 
(h = 0) BCS value A BCS = % e -t ^-i- 1 ^-^)/^^ 1 ' 2 
and a maximum at A = 0. For a higher Zeeman 
field A^cs/2 < h < AbcsI v2 ; the normal state at 
A = becomes a local minimum separated from the 
/i-independent global minimum at Abcs by a maxi- 
mum at A Sarma = Abcs y/2h/ A B cs - 1 0- For 
/i > Abcs/\/2 the minimum at A = lowers below 
that of the BCS state. For a fixed /i, this predicts a first- 
order SF-N transition at /i c (/x,<5) , with asymptotic form 
in the narrow FR limit given by 
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where ai, 2 = 8e- 2 /y/2, 1 (120 2/5 e" 8/5 , 4/5), for /i < 
(5/2 (/x 3> <5/2). The transition is accompanied by a jump 



in atom density from 71^(^1, 8) w |A/"(/i)/U + 2g 2 A 2 



down to nW(M,fc c ) = §c[0*+ ^c) 3/2 + (m~ M 3/2 ©(m 
/i c )] « n (s) (^i, <5) - pr/i 2 (l -7(m)/8), a jump in species 
imbalance from to to ~ 2J\f(fj,)h c , as well as other stan- 
dard thermodynamic singularities. 

In a more experimentally relevant ensemble of fixed 
total atom number n = —dEc/d/i, for h c \ = 
h c (n (s *>(n,8),5) < h < h c2 = h c (^ N \n,h),8) neither 
SF nor N states can satisfy the atom number constraint 
while remaining a ground state; ^- S - N ^> are SF and N 
chemical potentials at density n, Zeeman field h and de- 
tuning (5, obtained by solving n = n^ S ' N '(fx^ S ' N ') above. 
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for 8 > 2e F , 



For a narrow FR, 7 <C 1, we find 

lg[n-lc(6/2)^} 1/2 , for£«2e F , 
iAjr^e 16 w^ «f / , for£»2e F . 
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where Aj? = A^cs^, e F ) and h'- N ' > {8/2,n) is the solu- 
tion of n = n (N \d/2,h^ N \). Hence for h cl < h < h c2 
the gas phase separates ||| into SF and N rich regions 
in 1 — x{h, 8) and x{h, 8) proportions, determined by 
the atom number constraint xv> N > + (1 — x)rS s ^ = 
n. In above rS s \^, c ,8) > n > rS N \^, c ,8) are the 
SF and N state densities computed along the criti- 
cal chemical potential fj, c {h,5) determined by Eq. 
with limiting values fJ, c {h c i.2,8) = fj.( s > N '{n, S). The 
fraction of the N state admixture is then given by 
x{h,8) = (n^^ c> S)-n)[n^^ c ,S)~n( N \fi c ,S)y 1 
ranging between and 1 for h c \ < h < h c i spanning 
the coexistence region. 

A single-valued relation between the magnetization 
(species imbalance) m{h,8) — ^c[(^i c (h, 8) + h) 3 ^ 2 — 
(fjt c (h, 8) — h) 3 / 2 Q(ii c — h)~\ and Zeeman field h in the nor- 
mal paramagnetic state allows us to reexpress above pre- 
dictions in terms the species imbalance number M = xm, 
that is, the quantity (rather than h) that we anticipate to 
be kept fixed in atomic gas experiments. As illustrated 
in Fig^in a phase diagram expressed in terms of m = 8n 
the fully-paired SF state is confined to the detuning axis 
(m = 0) and the boundary between the coexistence re- 
gion and the N state is given by m C 2(n, 8) = m{h C 2, 8). 

We now turn to the negative detuning (BEC) regime. 
Although Eq, Eq. (J2J) and the phase diagram that follows 
from it can be accurately computed numerically (Fig^l, 
considerable insight can be gained by analytical analysis. 
This is particularly simple in the 7 — > limit in which 
E G ra {8-2 f i)\b\ 2 -{<ic/l5){h + n) 5 / 2 0{h + n),m{h, f i) = 
{2c/ 3) {h + ^) 3 ' 2 Q{h + n), and n = 2\b\ 2 + to. 

For h = and 8 < 0, this shows that the BCS su- 
pcrffuid ground-state smoothly crosses over to a BEC 
of closed-channel molecules, with a finite atom excita- 
tion gap yV 2 + g 2 |&| 2 ~ enforcing atom vacuum and 
the condensate density \b\ 2 w n/2 + 0(7). The gap 
equation dEo/db = 2{8 — 2y)b w then determines 
fx w (5/2, as in the crossover region, e^ 1 / 2 < 8 < 2e F , 
above. From the excitation spectrum E^ a it is clear that 
this ground state remains stable for < h < h m {8), 
with h m {8) = yJg 2 n/2+ {8/2) 2 w -5/2 determined by 
Eo,i{8, h m ) = 0. However, for h > h m , finite species im- 
balance to ps (2c/3)(ft — |(5|/2) 3 / 2 develops, depleting the 
SF condensate \b\ 2 = n/2 - {c/3){h - \S\/2) 3 / 2 . The re- 
sulting magnetized SFm is stable for h m < h < h C 2, with 
h c2 {8) = (3n/2c) 2 / 3 + \S\/2 = 2 2 / 3 e F - 6/2 determined 
by m(h C 2, 5/2) = n, giving a smooth extension into the 
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FIG. 2: Bogoliubov sound speed u in the BEC regime as a 
function of species imbalance m for S — — 2ef, vanishing at 
boundary of the SFm with the SF — N coexistence region. 

BEC regime of Eq. J7J, computed inside the BCS and 
crossover regimes. 

For a narrow FR Eq (A, fi, 5, h) can be accurately com- 
puted analytically giving 
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m{h')dh'. (8) 



Its minimization together with the atom number con- 
straint fixes fi » 8/2 + 0(7) and leads to the phase 
diagram in Fig. ^ We find[l5j that above expressions 
for h m {8) and h C 2{8) receive only small 0(7) correc- 
tion for 8 < — 7 1 / 2 e F . Hence in contrast to the BCS 
side, where the system undergoes phase separation for 
an arbitrary small m ^ 0, on the BEC side the tran- 
sition at h m {8) is into a uniform magnetized superfluid 
(SFm) that persists over a finite range of to and is a 
coherent superposition of a singlet molecular conden- 
sate and fully spin-polarized Fermi gas. The sequence 
SF^ SFm — > N of continuous transitions remains un- 
changed for 8 < 8 C ~ — 10.6e F . However, for a finite 7 
and 8 C » — 10.6e F < 8 < — 7 1 / 2 e F a secondary local (N 
state) minimum develops at A = leading to a lst-order 
SFm transition at h c i{5) » —0.658 < h C 2{8), pre- 
empting a continuous one at h C 2{8). For a fixed atom den- 
sity n and h c i{5) < h < h C 2{5) the gas phase-separates 
into coexisting SFm and N states. 

This h c i{8) boundary (equivalent to to c1 (<5) = 
m{h cl {6),6/2) w 0.029n|(5/e F | 3 / 2 ) is accurately (to 
C(7 2 )) located by the vanishing of the coefficient of A 4 
in Eq(A), proportional to the effective molecular scat- 

tering length a m {h,8) = effigy *WH)- We then 
predict 0| that the Bogoliubov sound velocity u{8,m) 
vanishes (to C(7), followed by a small jump to 0) at the 
lst-order SFm-N transition and exhibits a — m 1 / 3 cusp 
singularity at the SF-SFm boundary. The full expression 
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is illustrated in Fig. [5] and given by (with 6 = S/ep) 



i r- / / 2 5 / 3 
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with F(x) = 1 
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a; tan 



and uo = jjj^f the sound velocity at m = 0. 

We now turn to the FFLO state. Because Q ^ pair- 
ing is driven by the mismatch of the up | and J, Fermi 
surfaces, with Q »s fept — fcF|, it is clear that the FFLO 
state can only be stable at large positive detuning. Com- 
puting Eq for Q 7^ to leading order in Aq -C (us- 
ing dimensionless quantities Aq = Ag/ep, ft = m/ £ f 
h = h/ep, Q = Qy/Ji/kp, and £g 



Ec/ce^ 2 ), we find: 
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At fixed /i, for a given 5 and ft, the ground state is deter- 
mined by minimizing £g(Aq, Q) over Aq (the gap equa- 
tion) and Q (equivalent to vanishing of the ground state 
momentum). For 5 > 2ep we find a lst-order SF-FFLO 
(preempting SF-N) transition approximately at ft c (<5, fi), 
Eq. J5J. At fixed atom number, for ft > ft c the gas phase 
separates into coexisting SF and FFLO states, approx- 
imately bounded above by h C 2{8), computed for Q = 
above. At slightly higher field, /ifflo^), we find that 
the FFLO state undergoes a continuous transition (that 
on general grounds we expect to be driven lst-order by 
fluctuations) into the N state. Numerical solution of the 
gap, number and momentum equations yields hpFL,o{$) 
(and thus TOfflo via Eq. ©, plotted in Fig.l) that in- 
terpolates between 0.754Ai?((5) for large S (in agreement 
with FF [lOj) and h c i(6,n) for 5 — + <5*, with the crossing 
point <5* « 2ep. This microscopically calculated value 
of S* > contrasts with the conclusion of Ref. H the 
latter based on a purely qualitative discussion, that has 
little quantitative predictive power, e.g., in determining 
the precise location of phases. 

In free expansion experiments, the FFLO state, most 
easily observed with a trap having a typical size that is 
large compared to Q^ 1 , should exhibit a BEC peak (ob- 
served by its projection onto the molecular condensate 0) 
shifted by hQ,t/m a (t expansion time) corresponding to 
the finite momentum ftQ(d) (Fig^b)) of its condensate, 
and a (spontaneous) Bragg lattice of pea ks in the more- 
likely case of multiple-Q pairing [H 111 El E3 Thc 
anisotropy of the FFLO pairing should also be reflected in 
"noise" experiments 18] sensitive to angle- dependence of 
pairing correlations across the Fermi surface. Our predic- 
tions of gapless atomic excitations in the SFm and FFLO 
states, as well as the vanishing of the molecular scattering 



length a m and of the 0th sound velocity u at the SF m~ 
N phase boundary should be observable through Bragg 
spectroscopy and reflected in thermodynamics (e.g., heat 
capacity that is power-law in T). We also expect stan- 
dard thermodynamic anomalies across phase transitions 
in Fignja), and phase separation accompanied by den- 
sity discontinuity and local density variation with detun- 
ing and atom imbalance across the coexistence region. 
Finally, because a gas trapped in a smooth potential 
V(r) is well-characterized by a local chemical potential 
fi(r) = /j,— V(r) (the Thomas-Fermi approximation), our 
fixed /i analysis is directly experimentally relevant. For 
negative detuning and finite species imbalance we pre- 
dict SF state in the cloud's core of radius ro(6), with 
density discontinuity to the outer-shell N state, with 
ro(<5) determined by ^t(ro) = fi c (5,m) (see center inset 
of Figda)). We expect that this shell structure should 
be readily observable, particularly if different hyperfinc 
states and closed and open channels can be imaged inde- 
pendently Details of these experimental predictions will 
be presented in a forthcoming publication ■ 
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